r = 10;
F = eye(r);
B = eye(r);
A = eye(r);
E = F;
H = diag(logspace(-4, 4, r));
M = B * ones(size(B, 2), r);
X = mess_sylvester_sparse_dense(A, 'N', H, 'N', M, E, F);